Mb-MV*  33t 


TECHNICAL 

LIBRARY 


AD 


TECHNICAL  REPORT  ARLCB-TR-83032 


ELECTRICAL  TRANSPORT  IN  LOW  RESISTIVITY 

AMORPHOUS  ALLOYS 


L.  V.  MEISEL 
P.  J.  COTE 


SEPTEMBER  1983 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  CENTER 

LARGE  CALIBER  WEAPON  SYSTEMS  LABORATORY 
BENET  WEAPONS  LABORATORY 
WATERVLIET  IN.Y.  12189 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


DISCLAIMER 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  positior  unless  so  designated  by  other  authorized 
documents. 

The  use  of  trade  name(s)  and/or  manufacture^ s)  in  this  report  does  not 
constitute  an  official  indorsement  or  approval. 


DISPOSITION 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  retun  it  to 
the  originator. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Pete  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER 

ARLCB-TR-83032 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

ELECTRICAL  TRANSPORT  IN  LOW  RESISTIVITY 

AMORPHOUS  ALLOYS 

S.  TYPE  OF  REPORT  4  PERIOD  COVERED 

Final 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR(a) 

-L.  V.  Meisel  and  P.  J.  Cote 

8.  CONTRACT  OR  GRANT  NUMBERfa) 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Armament  Research  &  Development  Center 

Benet  Weapons  Laboratory,  DRSMC-LCB-TL 

Watervliet »  NY  12189 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 

AREA  4  WORK  UNIT  NUMBERS 

AM  CMS  N0.61H.02.Hb00.011 

PR0N  NO. 1A325B541A1A 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Armament  Research  &  Development  Center 

Large  Caliber  Weapon  Systems  Laboratory 

Dover*  NJ  07801 

12.  REPORT  DATE 

Sept.  1983 

13.  NUMBER  OF  PAGES 

33 

14.  MONITORING  AGENCY  NAME  &  ADDRESS^//  different  from  Controlling  Office) 

IS.  SECURITY  CLASS,  (of  this  report) 

Unclassified 

ISa.  DECLASSI  FI  CATION/ DOWN  GRADING 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thle  Report) 

Approved  for  Public  Release , Distribution  Ullimited 

17.  DISTRIBUTION  STATEMENT  (of  the  ebatract  entered  In  Block  20,  It  different  from  Report) 

is.  supplementary  notes  Presented  at  the  meeting  of  American  Physical  Society, 

March  1983.  Submitted  to  the  Journal ,  Physical  Review. 

19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  Identify  by  block  number) 

Electrical  Transport 

Amorphous  Alloys 

Diffraction  Model 

Saturation  Effects 

20,  ABSTRACT  fCarrtfnum  m  reverse  si rtm  ft  necreasary  and  Identify  by  block  number) 

Diffraction  model  calculations  incorporating  appropriate  scattering  matrix 
elements  and  phonon  ineffectiveness  effects  (saturation  effects)  yield  results 
which  are  consistent  with  the  observed  temperature  dependence  of  the  electrical 
resistivity  in  low  resistivity  (p  <  100  yftcm)  amorphous  alloys.  In  particular, 
remarkably  good  quantitative  agreement  with  available  detailed  resistivity 
measurements  in  a-MgZn  alloys  has  been  obtained  by  these  methods.  The  results 

( CONT 1 D  ON  REVERSE) 

DD  .UST,  1473  EDITION  OF  1  NOV  65  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dete  Entered ; 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfHTian  Dala  Entered) 


20.  ABSTRACT  (CONT'D) 

further  indicate  that  saturation  effects,  which  dominate  the  temperature 
dependence  of  the  high  resistivity  amorphous  metals,  are  important  even  for 
resistivities  as  low  as  50  yftcm. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE(T*7ien  Data  En tarad) 


TABLE  OF  CONTENTS 


Page 


ACKNOWLEDGEMENTS 

INTRODUCTION 

THEORY 

THE  SCATTERING  MATRIX  ELEMENT  AND  THE  MAGNITUDE  OF  p 
COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 
SUMMARY  AND  CONCLUSIONS 
REFERENCES 


TABLES 


I.  SUMMARY  OF  INPUT  FOR  RESISTIVITY  CALCULATIONS  IN  a-MgZn.  THE 
DEBYE  TEMPERATURE  0  WAS  DETERMINED  FROM  SPECIFIC  HEAT  MEASURE¬ 
MENTS,40  THE  POSITION  OF  THE  FIRST  PEAK  IN  THE  STRUCTURE 
FACTOR  kp  WAS  DETERMINED  BY  NEUTRON  DIFFRACTION, 39  AND  THE 
FERMI  WAVENUMBERS  kF  WERE  DEDUCED  FROM  ALL  EFFECT  MEASURE¬ 
MENTS.24  n  IS  THE  HARD  SPHERE  PACKING  FRACTION  IN  THE 
PERCUS-YEVICK  FORMULA.37 


II. 


A  SUMMARY  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS  IN  LOW 
RESISTIVITY  AMORPHOUS  ALLOYS.  THE  "BRILLOUIN  SCATTERING 
PART”  REFERS  TO  THE  PART  OF  THE  RESISTIVITY  ARISING  FROM 
THE  DELTA  FUNCTION  AT  K  =  0  IN  a(K).  THE  VALUES  OF  2kp 
AND  kD  AND  HENCE  2kF/kp  IN  THE  a-MgZn  ALLOYS  WERE 
EXPERIMENTALLY  DETERMINED;  THE  VALUES  OF  2kp/kp  GIVEN 
FOR  a-CuSn  ARE  NFE  VALUES  AND  ARE  CONSISTENT  WITH  THE 


CONCENTRATION  DEPENDENCE  OF  p  ;  THE  CONCENTRATION 
DEPENDENCE  OF  p  IN  a-AuSn  SUGGESTS  THAT  2kF/kp  IS 
~  10%  less  THAN  NFE  VALUES.  THE  COEFFICIENT  A  IS 
DERIVED  BY  p  «  p0(l+AT2)  FOR  5K  <  T  <  40K.  THE 
COEFFICIENT  B  IS  DEFINED  BY  p  «  Pi  -  BT3'2  FOR  %  < 

T  <  240K.  THE  COLUMN  HEADED  AVG  LISTS  AVERAGE  VALUES 
FOR  THE  ALLOYS  STUDIED  IN  REFERENCE  24. 


13 


14 


t 


i 


Page 


* 


LIST  OF  ILLUSTRATIONS 


1.  Weighted  absolute  squared  potentials  used  to  describe  transport 
in  a-MgZn. 

10 

2.  Electrical  resistivity  vs.  temperature  below  70K  in  a-MgZn  from 
Reference  24. 

19 

3.  The  relative  change  in  the  resistivity  vs.  T/ 0  computed  from  the 
diffraction  model  with  the  adjusted  phase  shift  expanded 
potential  appropriate  to  a-MgZn. 

20 

4.  The  3/2  power  region  to  the  right  of  the  resistivity  maximum  in 
a-MgZn. 

24 

% 


#' 


ii 


c 


ACKNOWLEDGEMENTS 


We  thank  Dr.  U.  Mizutani  and  Dr.  T.  Matsuda  for  their  helpful  comments 
and  for  providing  data  on  a-MgZn  prior  to  publication.  We  are  also  pleased  to 
acknowledge  the  assistance  of  Ellen  Fogarty  in  preparing  the  manuscript. 


iii 


INTRODUCTION 


Most  quantitative  theoretical  results  on  electrical  transport  in 
amorphous  metals  have  been  obtained  using  the  diffraction  model and  its 
extensions in  which  the  electron-phonon  interaction  matrix  element  is 
assumed  to  be  independent  of  the  electron  mean  free  path.  We  shall  refer  to 
such  calculations  in  this  report  as  applications  of  the  standard  diffraction 
model.  The  standard  diffraction  model  also  serves  as  the  basis  for  the 
analysis  of  transport  in  crystalline  metals,  yielding  Bloch— Gruneisen  theory^ 
for  normal  scattering  processes.  However,  the  standard  diffraction  model  is 
of  questionable  validity  for  most  amorphous  metals  whose  resistivities  are  of 
the  order  of  or  greater  than  150  pi2cm,  which  corresponds  to  electron  mean  free 
paths  of  the  order  of  ionic  spacings.  The  deviations  of  experimental  results 
from  the  predictions  of  the  standard  diffraction  model  in  high  resistivity 
metals  are  referred  to  as  "saturation  effects"  or  Mooij  phenomena. 1 3-17  Xn 
spite  of  the  presence  of  saturation  effects,  the  standard  diffraction  model 
gives  reasonable  values  for  the  magnitude  of  the  electrical  resistivity  and 
its  concentration  dependence  in  a  number  of  amorphous  alloys  ;6>18  in  fact,  the 
model  is  often  not  observed  to  fail  unless  close  attention  is  paid  to  details 
of  the  temperature  dependence  of  the  resistivity.  We  refer  to  treatments  of 
electrical  transport  in  which  the  electron-phonon  interaction  is  electron  mean 
free  path  dependent 1 6, 17  as  "saturated  cases"  of  the  diffraction  model. 


References  are  listed  at  the  end  of  this  report. 
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The  treatment  of  tfooij  phenomena  has  been  the  subject  of  considerable 
theoretical  study.  Some  investigators  have  treated  high  resistivity  metals  in 
the  context  of  the  diffraction  model  by  postulating  interband  tunneling 
channels  or ,  in  analogy  with  the  Pippard  theory  of  ultrasonic  attenuation  in 
metals ,20  an  electron  mean  free  path  dependent  electron— phonon  inter¬ 
action.16*17  Other  investigators21  have  approached  the  problem  by  extending 
theories  intended  for  transport  in  materials  whose  electrons  are  localized. 

The  standard  diffraction  model  is  expected  to  be  valid  when  the  electron 
mean  free  path  is  not  too  short.  Thus,  a  test  of  the  theory  in  low 
resistivity  (p  <  100  pflcra)  metallic  glasses  is  of  interest.  Until  recently, 
relatively  few  experiments  have  been  performed  on  such  alloys  since  they  are 
difficult  to  fabricate.  Notable  exceptions  were  the  results  in  vapor 
deposited  a-CuSn22  and  a-AuSn23  alloys  for  a  wide  range  of  compositions. 

These  alloys  exhibit  the  features  predicted  by  the  diffraction  model  including 
trends  in  p  and  the  temperature  coefficient  of  resistivity  (TCR)  with 
composition,  and  the  presence  of  a  maximum  in  p  vs .  T  in  alloys  with  a 
negative  room  temperature  TCR,  an  essential  feature  of  standard  diffraction 
model  predictions.6 

There  are  now  substantial  data  for  low  resistivity  amorphous  alloys. 
Matsuda  and  Mizutani  completed  a  thorough  study  of  electrical  transport  in 
a-MgZn  alloys2^  from  2  to  300K  with  Zn  concentrations  ranging  from  22.5  to  35 
percent  and  in  an  a— MgCu  alloy.23  Mizutani  and  Yoshida26  provided  data  on  a 
variety  of  Ag— Cu  based  alloys  from  77  to  300K.  Earlier ,  Hafner  et  al27  had 
found  good  agreement  between  experiment  and  the  theory  of  References  7-9  in  a 
preliminary  study  on  a-Mg7Zn3#  These  alloys  are  particularly  suited  for  tests 
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of  the  diffraction  model  because  of  their  low  resistivities  (~  50  yficm)  and 
the  existence  of  supporting  data  which  determine  the  relevant  parameters  in 
the  Ziman-Faber  theory  (e.g.,  Fermi  wave  number,  kp,  and  structure  factor  peak 
position,  kp) .  Moreover,  the  conduction  electron  states  in  a-MgZn  may  be 
assumed  to  be  almost  exclusively  of  s  and  p  character  in  contrast  to  the 
complex  situation  in  the  more  common  transition-metal  based  glassy  alloys. 

In  this  report  we  compare  the  a— MgZn  alloy  transport  data  with  results 
computed  using  the  diffraction  model  with  (i)  Heine-Aberenkov  pseudopotentials 
as  tabulated  by  Harrison28  used  for  the  scattering  matrix  element  in  the  Born 
approximation,  and  (ii)  a  phase  shift  expansion  of  the  scattering  matrix 
elements  as  developed  by  Evans  et  al.29  Saturation  effects  are  taken  into 
account  by  invoking  the  Pippard— Ziman  constraint  on  the  electron  phonon 
interaction. 1,6,16, 17  a.  brief  review  of  the  theory  is  given  in  the  section 
below.  The  scattering  potential  is  discussed  in  the  section  following  the 
Theory,  where  the  method  of  selecting  phase  shifts  appropriate  for  a— MgZn  is 
described  ;  and  the  problems  encountered  with  the  Born  approximation  are 
discussed.  In  the  section  on  Theoretical  and  Experimental  Results,  detailed 
comparison  is  made  between  the  data  on  a— MgZn  and  the  predicted  results  using 
the  diffraction  model  with  both  the  phase  shift  expansion  and  pseudopotential 
scattering  matrix  elements.  Implications  for  transport  in  general  low 
resistivity  alloys  are  also  described  in  that  section.  A  summary  and 
conclusions  are  given  in  the  last  section. 
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THEORY 


The  diffraction  model  (Ziman-Faber  theory^)  result  for  the  electrical 
resistivity  is 

12^o  ,1  K  K 


e2hVF2  '0  2kF  2kF 

where  is  the  atomic  volume,  Vp  the  Fermi  velocity,  kF  the  Fermi  wave 
vector,  K  the  scattering  vector,  h  is  Planck*  s  constant  divided  by  2tt  ,  e  is 
the  electron  charge,  and  the  resistivity  static  structure  factor  SP(K)  is 
defined  in  terms  of  the  Van  Hove  dynamical  structure  factor^  S(K,w)  as 


/  d( - )( - )  SP(K)|u(K)|2 


(1) 


SP(K)  =  J  dwxn(x)  S(K,  cj)  (2) 

— oo 

where  x  -  hoo/kg,  n(x)  =  (ex-l)*“1  ,  kg  is  Boltzmann’s  constant,  and  T  is  the 
absolute  temperature.  The  scattering  matrix  element  u(K)  is  approximated  by  a 
pseudopotential  in  Born  approximation  or  for  strong  scattering  (i.e.,  large 
phase  shifts) ,  it  can  be  expressed  in  terms  of  phase  shifts2^  as 

27rh3  x-  inn(E^) 

u(K)  - - - l( 2Jl+l)sin  n*(EF)e  *  F  P*(cos  6)  (3) 

m(2mEF)1/2^0  i 

where  the  phase  shift  n^(EF)  for  angular  momentum  quantum  number  %  is 
evaluated  at  the  Fermi  energy  and  ra  is  the  electron  mass.  The  matrix  element 
in  the  form  of  Eq.  (93)  includes  single  site  multiple  scattering  and  is  called 
the  t -matrix.  We  shall  discuss  transport  in  a-MgZn  alloys  using  both  of  these 
formulations.  The  general  case  is  described  in  terms  of  the  t -matrix. 

The  resistivity  static  structure  factor  in  an  amorphous  Debye  solid  may 
be  expanded  in  the  form^”^ 

SP(K)  =  SQ  P ( K)  +  SXP(K)  +  S2p(K)  +  ...  (4) 
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where  (K)  is  an  n-phonon  term.  The  elastic  scattering  term  (no  phonons) 
is 

SQ  P  (  K)  =  a(K)e“2W(K)  (5) 

where  e”2w(K)  is  the  Debye-Waller  factor  and  the  geometrical  structure  factor 


a(K)  is 


a(K) 


1 

N 


I  exp [ IK* (m-n) ] 
m  ,n 


(6) 


with  m,n  averaged  ionic  positions.  The  one  phonon  term,  allowing  for  Pippard- 
Ziman  phonon  ineffectiveness  as  described  in  References  16  and  17,  is 


✓  T  ,  1  q  q  o  ,  d  ftq  ,  —  — . 

SiP(K)  =  a(K)e~2W(K)  -  j  d(— )  (— )  n(x)(n(x)+l)  / - a(  |K+q|)F(qA)  (7) 

0  0  qD  qD  4it 

where  ct(K)  =  3(hK)2/Mk.B0  where  M  is  the  averaged  ionic  mass,  qp  is  the  Debye 
wave  number,  0  is  the  Debye  temperature,  x  =  (e/TXq/qp)  for  a  Debye  solid,  A 
is  the  electron  mean  free  path,  and  F(qA)  describes  the  reduction  in 
scattering  effectiveness  which  occurs  for  small  qA.  The  calculations 
presented  here  assume  that  F(qA)  can  be  represented  by  the  form  suggested  by 
Pippard^O  in  his  study  of  ultrasonic  attenuation,  and 

2  ytan-1y  3 

F(y)  =  -  [ - -  -]  .  (Pippard  saturation)  (7a) 

y-tan^y  Y 

We  refer  to  this  expression  as  the  Pippard  function#  If  saturation  effects 
can  be  ignored 

F(y)  =  1#  (No  saturation)  (7b) 

We  refer  to  this  form  of  the  theory  as  the  standard  diffraction  model. 

Equation  (7b)  is  also  clearly  the  long  electron  mean  free  path  limit  (or  low 
resistivity  limit)  of  the  Pippard  function.  One  can  also  represent  saturation 
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effects  with  a  "sharp  cutoff"  form 

0  y  <  yc 

F(y)  *  (sharp  cutoff)  (7c) 

i  y  >  yc 

where  the  cutoff  value  yc  is  of  order  qpa  where  a  is  the  mean  ionic  spacing. 
Note  that  the  placement  of  F(qA)  in  the  one  phonon  part  of  the  (generalized) 
resistivity  static  structure  factor  results  from  a  decomposition  of  the 
absolute  square  of  a  phonon  wave  number  and  electron  mean  free  path  dependent 
electron-phonon  matrix  element  into  a  simple  product  according  to  the 
prescription 

|u(K,qA)|2  =  |u(K) | 2F(qA) 

where  u(K)  is  the  ordinary  scattering  matrix  element  which  appears,  for 
example,  in  the  elastic  scattering  term. 

The  multiphonon  series  is  approximated  in  various  ways. 30  For  the  range 
of  temperatures  of  interest  here  (T  <  0)  the  particular  approximation  made  is 
not  important  and  the  results  are  given  in  Sham— Ziman  approximation, 30  i.e., 
we  assume  that  the  effect  of  the  multiphonon  series  is  to  cancel  the  Debye- 
Waller  factor  in  the  one— phonon  term.  We  shall  refer  to  that  modified  terra  as 
the  inelastic  scattering  term  or  simply  the  phonon  scattering  term. 

The  Debye-Waller  exponent  is  given  (for  Debye  solid)  by 

1  q  q  1 

2W(K)  =  a(K)J  d( — )  ( — )  (n(x)+  -)  (8a) 

0  qD  QD  2 

T  2  0/T  1 

=  ct(K)(-)  /  dx  x(n(x)+  -)  (8b) 

0  0  2 
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Note  that  Eqs.  (1)  through  (8)  have  been  written  in  a  fora  appropriate  to 
a  pure  substance.  The  product  SP(K)|u(K)|2  should  be  replaced  by  a  sum  of 
concentration  dependent  terras  involving  individual  constituent  scattering 
matrix  elements  and  partial  structure  factors  in  alloy  systems.  Nevertheless, 
in  the  a-MgZn  alloys  we  treat  the  material  with  a  single  effective  scattering 
potential  in  the  spirit  of  the  substitutional  model  of  Faber  and  Ziman^  and  a 
Percus-Yevick  formal  for  the  geometrical  structure  factor.  This  is  not  quite 
right  when  there  is  short  range  order  and  2kp  is  situated  differently  with 
respect  to  the  peak  positions  of  the  various  partial  structure  factors.  How¬ 
ever,  a  broad  range  of  scattering  vectors  contribute  to  the  resistivity  in  a- 
MgZn,  the  constituents  have  very  similar  atomic  structure,  and  the  computa¬ 
tions  by  von  Heimendal^  suggest  that  the  three  partial  structure  factors  are 
very  similar;  we  thus  assume  that  the  conditions  for  the  substitutional  model 
obtain.  The  application  of  the  substitutional  model  to  the  other  low  resis¬ 
tivity  amorphous  alloys  considered  here  may  not  be  as  good  an  approximation. 

THE  SCATTERING  MATRIX  ELEMENT  AND  THE  MAGNITUDE  OF  p 

The  atomic  structure  of  Mg  consists  of  filled  shells  and  3s2  electrons 
and  similarly  in  Zn,  filled  shells,  and  4s2  electrons.  It  follows  from 
straightforward  considerations  of  the  atomic  structure  of  these  simple 
divalent  metals  that  the  appropriate  starting  structure  in  the  metallic  solid 
will  consist  of  a  nearly  filled  s  band  with  some  p  band  occupation  resulting 
in  a  metallic  sp  band  primarily  of  s-character  with  a  small  admixture  of 
p-character.  (Note  that  the  energy  required  to  promote  an  electron  to  a 
d-level  in  Mg  or  Zn  is  prohibitively  large  in  constrast  to  the  situation  in 
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Ca,  Sr,  or  Ba.)  Thus,  if  we  construct  a  scattering  matrix  element  as  given  in 
Eq.  (3),  we  expect  the  s-wave  phase  shift  nQ(Ep)  to  be  slightly  less  than  tt  , 
the  p  phase  shift  rii(Ep)  to  be  small,  and  other  phase  shifts  to  be  negligible. 
Such  a  scattering  matrix  element  will  be  drastically  different  in  form  from  a 
Born  approximation  pseudopotential  matrix  element.  We  have  constructed  such  a 
matrix  element  for  a-MgZn  by  adjusting  no(EF)  to  Sive  the  observed  magnitude 
of  p  with  m(Ep)  constrained  to  satisfy  the  Friedel  sum  rule. 33  (Since  there 
are  only  &  =  0  and  1  phase  shifts  to  consider  and  the  electron  per  atom  ratio, 
z  =  2,  one  has  to  satisfy  uq  +  3ni  =  tt.)  The  resulting  values  of  tiq(ef)  and 
rii(Ep)  are  2.87  and  0.09  respectively;  we  use  these  values  in  our  subsequent 
investigations  of  the  temperature  dependences. 

Dunsworth34  has  deduced  phase  shifts  for  Zn  in  8  brasses  by  adjusting 
APW-form  matrix  elements  to  fit  experimentally  determined  Fermi  surface 
features.  His  results  for  Zn  in  crystalline  8  brasses  where  Z  =  1.5  electrons 
per  atom  are  n0  =  3.114,  ni  =  0.289,  and  H2  =  0.001  at  the  free  electron  Fermi 
energy.  This  lends  support  to  our  general  expectation  that  the  divalent 
metals  and  alloys  will  exhibit  large  l  =  0  phase  shifts  (near  tt)  and  small  &  = 
1  phase  shifts. 

It  should  be  stressed  that  although  the  magnitude  of  p  is  quite  sensitive 
to  changes  in  the  phase  shifts,  the  T  dependence  of  p(T)/p(0)  is  relatively 
insensitive  to  specific  phase  shift  values.  For  example,  a  matrix  element 
constructed  for  nQ  -  2.00  and  =  0.38  yields  a  factor  of  two  increase  in  p 
but  only  changes  the  temperature  dependence  of  p(T)/p(0)  by  about  5  percent 
over  the  range  from  0°K  to  0  for  the  values  of  2kp/kp  studied  here.  In  fact, 
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TT 

for  -  <  nQ  <  tt  and  ni  adjusted  to  satisfy  the  sum  rule,  the  scattering  matrix 

2 

elements  are  similar  in  form  to  that  shown  in  Figure  1  for  ni  =  2.87  and  give 
essentially  equivalent  (i.e.  ,  within  about  20  percent)  temperature  dependences 
for  p(T)/p(0) . 

Another  point  which  should  be  considered  is  the  effect  on  p(T)  of 
breakdown  of  the  Freidel  sum  rule  constraint.  Although  various  forms  of  the 
Friedel  sum  rule^S,  36  are  invoked  to  constrain  phase  shifts  evaluated  at  the 
Fermi  energy,  it  should  be  noted  that  calculated  phase  shifts  often  fail  to 
satisfy  the  sum  rule;  for  example,  in  Reference  36  the  phase  shifts  computed 
for  Ca  and  Sr  yield  Friedel  sums  differing  by  13  and  3  percent  respectively 
from  the  sum  rule.  Again  the  form  of  the  t-matrix,  (which  determines  the 
temperature  dependence  of  p(T)/p(0))  is  preserved  even  for  potentials  which  do 

TT 

not  satisfy  the  sum  rule  as  long  as  -  <  ri0  <  tt  and  is  not  too  large. 

The  point  of  this  discussion  is  that  although  the  phase  shifts  deduced 
for  a-MgZn  are  not  uniquely  determined  by  the  transport  data,  the  temperature 
dependence  of  the  resistivity  of  all  the  low  resistivity  amorphous  alloys 
studied  to  date  can  be  discussed  in  terras  of  an  approximate  t-matrix  of  the 
form  shown  in  Figure  1.  Also  one  might  expect  the  best  phase  shifts  for 
a-MgZn  to  be  given  within  about  the  deviations  from  the  sum  rule  seen  in  Ca  or 
Sr  • 

If  one  uses  the  Born  approximation  with  Heine-Aberenkov  pseudopotentials 
as  tabulated  in  Reference  28  for  Hg  and  Zn,  approximates  SP(K)  by  S0(K) ,  and 
takes  a  value  of  kp  appropriate  to  the  a-MgZn  alloys,  then  one  obtains  20  and 
^  29  uftcm  respectively  for  the  I^g  and  Zn  pseudopotentials.  The  discrepancy 
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Figure  1.  Weighted  absolute  squared  potentials  used  to  describe  Transport  in  a-MgZn. 


between  either  of  these  values  and  the  measured  resistivities  of  the  a-MgZn 
alloys  is  of  the  order  found  when  these  methods  have  been  applied  to  transport 
in  liquid  metals • 37  However,  the  form  of  the  Born  approximation  matrix 
element  is  drastically  different  from  that  of  the  adjusted  phase  shift  matrix 
element  as  may  be  seen  in  Figure  1.  The  cancellation  near  2kp ,  which  is 
characteristic  of  pseudopotentials,  does  not  occur  in  the  phase  shift  expanded 
matrix  element. 

Considering  only  the  theoretical  results  for  the  magnitude  of  p,  either 
approach  can  yield  reasonable  agreement  with  the  data  in  a-MgZn  alloys.  How¬ 
ever,  we  shall  see  in  the  section  on  Results  that  the  phase  shift  matrix 
element  leads  to  a  better  approximation  to  the  observed  temperature  dependence 
in  a-MgZn  than  the  Born  approximation  results. 

Moreover,  there  are  well  known  examples  of  difficulties  with  Born 
approximation  transport  calculations  in  related  systems.  For  example,  the 
electrical  resistivity  computed  with  Heine-Aberenkov  pseudopotentials  in  Born 
approximation  and  with  Percus— Yevick  hard  sphere  structure  factors  is  in  poor 
agreement  with  experiment  in  liquid  Ca,  Ba ,  and  Sr.  On  the  other  hand,  Ratti 
and  Evans33  obtained  good  agreement  with  the  experimental  data  in  liquid  Ca , 
Ba,  and  Sr  using  a  phase  shift  expanded  form  for  the  scattering  matrix  element 
(Eq.  (3))  with  phase  shifts  computed  from  muffin  tin  potentials.  (The  only 
significant  phase  shifts  in  these  alkaline  earth  liquids  were  an  s  phase  shift 
slightly  less  than  tt  and  a  small  d  phase  shift.) 

Furthermore,  there  are  even  indications  of  difficulties  with  Born 
approximation  for  the  treatment  of  transport  in  monovalent  liquid  metals. 

Young  et  al33  concluded  that  Born  approximation  was  inadequate  in  these  liquid 
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metals  and  achieved  improved  agreement  with  electrical  resistivity  and  thermo- 
power  in  terms  of  phase  shift  expanded  matrix  elements  (Eq.  (3)).  In  the 
light  of  these  findings  for  liquid  Ca,  Ba ,  and  Sr  and  for  the  liquid  mono¬ 
valent  metals,  it  is  perhaps  not  surprising  that  we  find  that  the  temperature 
dependence  of  the  electrical  resistivity  in  a— MgZn  and  the  other  low 
resistivity  amorphous  alloys  is  better  described  by  a  scattering  matrix 
element  of  the  form  of  Eq.  (3)  (t-matrix)  than  by  pseudopotentials  in  Born 
approximation. 

COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS 

The  parameters  used  in  the  computations  are  listed  in  Table  I.  They  were 
based  upon  the  following:  x-ray  and  neutron  diffraction^  yield  values  of  kp 
in  a-MgyZn3  of  2.7  and  2.6  A-1  respectively.  Hall  effect  measurements  by 
Matsuda  and  Mizutani2^  give  2.77  A-1  <  2kp  <  2.98  A-1  for  the  range  of 
compositions  of  a-MgZn  studied.  The  mean  ionic  mass  is  used  for  M  and  in 
a-MgZn  we  take  qp  =  kp  (since  z  =  2) .  Specific  heat  measurements^®  in 
a-Mg7Zn3  yield  6  =  29 5K.  We  have  also  assumed  a  Debye  phonon  spectrum  and  an 
effective  (geometric)  structure  factor  of  Percus— Yevick  form  with  packing 
fraction  n  =  0.525  (which  is  representative  of  structure  factors  found  in 
amorphous  alloys) • 

In  the  remainder  of  this  report,  when  we  refer  to  the  saturated  case,  we 
are  quoting  results  computed  for  Pippard  saturation  with  qpA  =  11.7  in  Eq. 
(7a),  which  produces  a  25  percent  reduction  in  the  inelastic  scattering  part 
of  the  resistivity  at  T  «  0.  (This  is  equivalent  to  assuming  the  saturation 
resistivity  p*  «  200  p^cm  in  the  treatment  of  Reference  17.)  We  emphasize 
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here  that  inelastic  scattering  contributes  less  than  five  percent  to  the  total 
p  so  that  the  effect  of  saturation  in  the  adjusted  phase  shifts  is  negligible. 
The  unsaturated  case  results  are  obtained  from  the  standard  diffraction  model, 
i.e.,  using  Eq.  (7b)  or  letting  A  *  00  in  Eq.  (7a). 

TABLE  I.  SUMMARY  OF  INPUT  FOR  RESISTIVITY  CALCULATIONS  IN  a-MgZn.  THE  DEBYE 
TEMPERATURE  0  WAS  DETERMINED  FROM  SPECIFIC  HEAT  MEASUREMENTS,40  THE 
POSITION  OF  THE  FIRST  PEAK  IN  THE  STRUCTURE  FACTOR  kp  WAS  DETERMINED 
BY  NEUTRON  DIFFRACTION, 39  AND  THE  FERMI  WAVENUMBERS  kF  WERE  DEDUCED 
FROM  HALL  EFFECT  MEASUREMENTS. 24  n  IS  THE  HARD  SPHERE  PACKING 
FRACTION  IN  THE  PERCUR-YEVICK  FORMULA.37 


Input 

Descriptive 

Parameters 

Phonon  Spectrum 

Debye  Form 

0  =  29  5K 
qD  =  kF 

Geometrical 

Percus-Yevick  Hard 

n  =  0.525 
kp  =  2.6  A-1 

Structure  Factor 

Sphere  Form 

Fermi  Wavenumbers 

For  Zn  Concentration  i 

For  Zn  Concentration  i 

0.225,  2kF  =  2.77  A-1 

0.35,  2kF  =  2.98  A-1 

Table  II  summarizes  many  of  the  principal  experimental  results  regarding 
electrical  transport  for  T  <  0  in  a-MgZn  and  also  in  a-AuSn  and  a-CuSn  alloys 
along  with  corresponding  theoretical  results  based  upon  Mg  and  Zn 
pseudopotentials  and  the  adjusted  phase  shift  matrix  elements  for  the 
saturated  and  unsaturated  cases.  Generally,  the  adjusted  phase  shift  matrix 
element  in  the  unsaturated  case  yields  better  agreement  with  the  detailed 
experimental  results  than  the  Born  approximation  results  and  the  phase  shift 
results  including  saturation  yield  the  best  agreement  with  experiment  in  the 
a-MgZn  alloys.  Furthermore,  there  is  excellent  qualitative  agreement  between 
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the  phase  shift  results  and  the  amorphous  noble  metal  based  data.  These 
results  strongly  suggest  that  saturation  effects  or  Mooij  phenomena  are 
readily  observable  in  amorphous  alloys  with  resistivity  as  low  as  50  yftcra. 

Let  us  now  turn  to  a  brief  description  of  our  results  for  the  various 
temperature  regions  studied  in  detail  by  Matsuda  and  Mizutani^  in  a-MgZn 
and  the  overall  temperature  dependences  predicted  for  the  other  known  low 
resistivity  amorphous  alloys, 

A,  The  room  temperature  (T  »  0)  TCR. 

The  high  temperature  (T  >  0)  form  for  the  (effective)  resistivity  static 
structure  factor  to  first  order  in  a(K)  is 

SP(K)  2  a(K)e-a(K)M  +  a(K)  (T/  0)  [ AP( K)  (  l-y)-a(  K)e“a(K) /4]  (9) 

where  y  =  0  in  the  standard  diffraction  model  (i,e,  no  saturation)  and  for 
Pippard  or  sharp  cutoff  saturation  y  53  p/p*,  Generally,  p*  200  y^cm  in  high 
resistivity  systems  14-17  ancj  so  we  have  taken  y  *  1/4  in  a-MgZn,  This  was  the 
basis  for  our  choice  of  q^A  =  11,7  in  Pippard  saturation.  The  high 
temperature  limiting  form  of  the  averaged  resistivity  structure  factor  AP(K) , 
as  defined  in  References  6-9  is  given  by 

AP(K)  =  /  d(— )  / - a(  | K+q  |  )  (10) 

0  qp  4tt 

The  sign  of  the  contribution  to  the  TCR  from  scattering  vector  K  is  determined 
by  AP(K)(l-y)  -  a( K)e“a(K) ,  In  the  extreme  backscattering  case  (as  is  often 
assumed  for  TM  based  amorphous  alloys)  one  has  p  a  SP(2kp)  and  the  TCR  is 
negative  if  a(  2kp)e”a(^^) >  AP(  2kp)  ( 1-y)  ;  in  particular,  for  no  saturation 
and  Percus-Yevick  structure  factors  with  packing  fraction  0,525  (appropriate 
for  a  large  class  of  amorphous  metals)  negative  TCR  would  be  predicted  for 
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2kF 

0. 9  <  -  <  1.1.  On  the  other  hand,  when  a  relatively  broad  range  of  scatter- 

kF 

ing  vectors  yield  significant  contributions  to  p  and  the  TCR,  one  must  employ 
Eq.  (1)  to  determine  p(T),  and  this  simple  criterion  for  the  occurrence  of 
negative  TCR  is  invalid.  Saturation  effects  will  also  invalidate  this 
criterion;  the  range  of  2kp/kp  for  negative  TCR  is  increased  when  y  >  0,  and 
negative  TCR  will  generally  occur  for  all  2kp/kp  values  when  y  >  0. 5. 

The  computed  TCR's  in  a-MgZn  are  listed  in  Table  II  for  the  various 
conditions  considered.  The  "Brillouin  scattering’*  contribution,  6,  10  which 
arises  from  the  delta  function  at  K  =  0  in  a(K)  and  gives  the  normal 
scattering  contribution  in  the  crystalline  case,  is  listed  separately  in  Table 
II  to  indicate  the  relative  importance  of  this  term.  (We  have  generally 
assumed  that  Brillouin  scattering  contributions  are  negligible  in  TM  based 
amorphous  alloys  since  backseat tering  is  expected  to  dominate.) 

A  number  of  interesting  features  are  indicated: 

(i)  The  range  of  2kp/kp  for  which  negative  TCRfs  are  predicted  with  the 
phase  shift  approximation  in  Eq.  (1)  is  shifted  toward  considerably  higher 
values  than  those  given  in  the  extreme  backscattering  case.  The  Mg  pseudo¬ 
potential  results  suggest  a  similar  but  smaller  shift.  The  range  of  2kp/kp 
for  negative  TCR  in  saturated  and  unsaturated  cases  (excluding  the  Zn  pseudo¬ 
potential  case)  are  consistent  with  the  data  of  Reference  24.  Negative  TCR 
values  were  obtained  for  0.96  <  2kp/kp  <  1.24  in  the  unsaturated  case;  the 
range,  when  saturation  effects  (qj)A  *  11.7)  are  included  is  0.94  <  2kp/kp  < 
1.29. 
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(ii)  The  computed  Brillouin  scattering  contribution  to  the  TCR  is 
negligible  in  the  phase  shift  formulation,  but  is  significant  in  the  Born 
approximation  case.  For  the  Zn  pseudopotential,  negative  TCR' 8  are  completely 
eliminated  by  the  Brillouin  contribution.  Frobose  and  Jackie10  had 
encountered  similar  difficulties  with  their  Born  approximation  calculations  on 
a-CuSn. 

( iii)  The  Born  approximation  results  for  the  TCR  are  in  poor  agreement 
with  experiment#  The  phase  shift  results  are  about  half  as  large  as  the 
observed  TCR. 

B.  The  low  temperature  T2  region,  minima,  and  maxima  in  p(T) . 

The  low  temperature  limiting  form  of  the  resistivity  static  structure 

factor, neglecting  saturation  is 

SP(K)  =  a(K)  +  (tt2/6)  a( K) a(K) (T/ 0) 2  OD 

Thus,  the  standard  diffraction  model  yields  a  positive  quadratic  temperature 
dependence  for  the  resistivity  near  0°K,  independent  of  the  sign  of  the  TCR  at 
high  temperatures.  Consequently,  a  general  feature  of  the  standard 
diffraction  model  is  that  if  the  room  temperature  TCR  is  negative,  the 
resistivity  will  exhibit  small  maxima  in  p(T) .  (This  simple  result  is  often 
ignored ,  leading  to  improper  conclusions  regarding  agreement  between  theory 

rs  / 

and  experiment  in  high  resistivity  amorphous  alloys.)  Matsuda  and  Mizutaniz 
observe  a  quadratic  p(T)  in  a- MgZn  below  30K,  which  they  attribute  to  the  low 
temperature  form  of  the  theory  based  upon  Eqs.  (1)  and  (11).  However, 
detailed  calculations  indicate  that  significant  deviations  from  Eq.  (11)  occur 
above  about  5K  for  K  near  kp  so  this  interpretation  is  questionable  ;  in 
particular,  the  low  temperature  limit  of  the  T2  coefficient  will  generally  not 
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be  observed  above  ~  5K.  Nevertheless,  the  standard  diffraction  model  results 
for  a-MgZn  are  consistent  with  a  quadratic  p(T)  in  the  range  from  5  to  30K. 
However,  the  computed  coefficient  of  the  quadratic  term,  given  in  Table  II,  is 
smaller  than  the  low  temperature  limit  (which  is  about  2xlCT6K”2)  in  the 
unsaturated  case  but  is  still  larger  than  the  measured  coefficient  by  a  factor 
of  two.  When  saturation  is  included,  excellent  agreement  with  the  measured 
coefficient  is  obtained. 

Figure  2  shows  the  p  vs,  T  curves  obtained  by  Matsuda  and  Mizutani^  for 
a-MgZn  alloys  between  2  and  7 OK,  Small  maxima  in  the  resistivity  (p(T^)  - 
p(0°K)  »  10“^  p(0°K))  are  observed  for  «  50K  where  is  the  temperature  at 
the  resistivity  maxima.  Figure  3(a)  shows  phase  shift  based  standard 
diffraction  model  results.  The  2kp/kp  -1,1  case  is  appropriate  for  a-MgZn 
and  is  in  qualitative  agreement  with  the  experimental  results  shown  in  Figure 
2.  However,  the  computed  maximum  occurs  at  T^  *  0,3  0,  which  corresponds  to 
T^  »  9 OK  and  the  computed  maximum  is  substantially  larger  than  observed.  When 
saturation  is  included,  excellent  agreement  (see  Table  II  and  Figure  3(b))  in 
the  position  and  size  of  the  maximum  is  obtained.  The  Born  approximation 
results  at  these  temperatures  are  in  very  poor  agreement  with  experiment. 
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Figure  3.  The  relative  change  in  the  resistivity  vs.  T/9  computed  from  the  diffraction 
model  with  the  adjusted  phase  shift  expanded  potential  appropriate  to  a-MgZn. 


The  data  of  Figure  2  exhibit  another  interesting  feature,  viz.  small 
minima  at  about  10K.  This  feature  is  not  consistent  with  standard  diffraction 
model  predictions.  However,  as  may  be  seen  in  Figure  3(b)  or  from  the  results 
listed  in  Table  II,  the  theory  including  saturation  is  again  in  excellent 
agreement  with  the  data.  This  is  a  particularly  interesting  result  since  to 
observe  such  effects,  qpA  has  to  be  large  enough  to  produce  an  observable 
minimum  but  not  so  large  as  to  produce  monatonic  decreasing  p(T)  ;  also  the 
alloys  under  investigation  must  be  free  of  effects  associated  with  magnetic 
ions  or  transitions  to  superconductivity.  (Essentially  the  same  results, 
although  with  deeper  minima,  are  obtained  in  the  case  of  sharp  cutoff 
saturation  corresponding  to  y  =  1/4.)  It  is  possible,  of  course,  that  this 

agreement  is  fortuitous. 

C.  Other  low  resistivity  amorphous  alloys. 

Other  low  resistivity  (p  <  100  pHcm)  amorphous  alloys  with  2kp/kp  >  1 
that  we  are  aware  of  and  for  which  the  temperature  dependence  of  the 
electrical  resistivity  has  been  determined,  are  a-CuSn,22  a-AuSn,23  a-AuIn,23 
a-MgCu,25,  a-AgCuAl,26  a-AgCuMg,26  and  a-AgCuGe.26  These  alloys  have  been 
studied  over  wide  ranges  of  temperature  and  composition  which  correspond  to 
extensive  ranges  in  2kp/kp  ;  for  example,  for  a- CuSn  and  a- AuSn,  this  range 
varies  from  about  1.0  to  about  1.3.  Figure  3  shows  theoretical  results  for 
similar  ranges  in  2kF/kp  in  the  saturated  and  unsaturated  cases  and  parameters 
appropriate  to  the  a- MgZn  alloys. 

Before  the  results  exhibited  in  Figure  3  are  compared  with  measurements 
in  the  other  alloys,  we  make  the  following  comments:  (i)  The  magnitude  of  the 
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variations  in  the  electrical  resistivity  with  T  in  the  range  from  0°K  to  0  is 
governed  essentially  by  M0.  Hence,  although  the  effective  masses  are  larger, 
the  expected  lower  Debye  temperatures  for  these  noble  metal  based  alloys 
relative  to  a-MgZn,  can  account  for  the  fact  that  the  temperature  dependent 
effects  in  the  resistivity  are  of  the  same  magnitude.  (ii)  The  range  of 
resistivities  in  the  other  alloys  is  considerably  larger  than  in  the  a— MgZn 
alloys.  Thus,  we  cannot  expect  qpA  =  11.7  to  be  appropriate  for  all  these 
alloys.  This  effect  can  be  significant;  e.g.,  for  p  ®  100  gftcra,  saturation 
effects  could  completely  eliminate  the  small  maximum  in  p(T) ,  for  2kp  ra  kp , 
yielding  a  monotonic  decreasing  function.  (In  fact,  this  was  observed  in 
Reference  23  for  the  a-AuSn  alloy  with  p  »  100  pi2cm.)  (iii)  There  could  be 
appreciable  short  range  order  effects  produced  by  the  differences  in  effective 
potential  and  application  of  the  "substitutional  model"  might  not  be 
appropriate  for  such  different  ionic  constituents  as  occur  in  these  alloys. 

In  spite  of  these  difficulties,  all  the  low  resistivity  alloy  systems 
exhibit  the  general  features  shown  in  Figure  3  including  resistivity  maxima  in 
negative  TCR  cases.  This  strongly  suggests  that  matrix  elements  of  the  form 
given  in  the  phase  shift  expansion  rather  than  Born  approximation  pseudo¬ 
potential  is  appropriate  in  these  alloys.  The  temperature  dependences  of  the 
various  composition  (i.e.,  various  2kp/kp)  a— CuSn  alloys  are  in  remarkable 
agreement  with  the  saturated  case  curves  in  Figure  3(b)  if  we  use  free  elec¬ 
tron  theory  to  compute  2kp/kp  for  each  composition.  Similar  consistency  is 
obtained  for  the  a-AuSn  data  if  we  deduce  2kp/kp  from  the  p  and  TCR  vs.  compo¬ 
sition  data  shown  in  Reference  23  (which  leads  to  2kp  values  about  10  percent 
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smaller  than  free  electron  values).  The  a-MgCu25  data  are  essentially 
identical  to  those  of  a— MgZn.  The  AgCu  based  alloys2^  offer  a  particularly 
clear  illustration  of  the  shift  to  higher  values  for  the  2kp/kp  range  of 
negative  TCR  values. 

D.  The  (T-Tm)3/2  temperature  dependence. 

Matsuda  and  Mizutani24 > 25  discovered  that  over  an  extensive  range  of 
temperature  to  the  right  of  T^,  the  resistivity  in  the  a-MgZn  and  a-MgCu 
alloys  has  the  form 

p(T)  =  Pi  -  B(T-Tm)3/2  (12) 

The  phase  shift  calculation,  including  saturation,  fits  this  equation,  and 
thus  the  experimental  data,  surprisingly  well  both  with  regard  to  the  value  of 
B  (Table  II)  and  the  linearity  of  the  p  vs.  (T-TM)3/2  as  seen  in  Figure  4. 

SUMMARY  AND  CONCLUSIONS 

Electrical  resistivity  of  low  resistivity  (p  <  100  pOcm)  amorphous  metals 
has  been  studied  in  the  context  of  the  Ziman-Faber  diffraction  model  which  has 
been  generalized  to  account  for  saturation  effects  by  incorporation  of  the 
"Pippard  function"20  which  describes  the  reduction  of  the  electron-phonon 
interaction  at  small  qA  where  q  is  the  phonon  wave  number  and  A  the  electron 
mean  free  path.  The  standard  diffraction  model  formulae  are  obtained  in  the 
limit  that  A  goes  to  infinity.  (Essentially  equivalent  results  were  obtained 
using  a  "sharp  cutoff"  form  to  describe  saturation  effects.) 
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Figure  4.  the  3/2  power  region  to  the  right  of  the  resistivity  maximum  in  a— MgZn 


The  specific  computational  results  were  obtained  for  a-MgZn  alloys  which 
comprise  an  ideal  system  to  test  the  diffraction  model  because  of  their 
relatively  simple  electronic  structure,  low  resistivity,  and  the  existence  of 
detailed  experimental  data  covering  temperatures  from  2  to  300K  in  a  series  of 
well  characterized  alloys.  Moreover,  the  similarity  of  the  atomic  structures 
and  ionic  radii  of  and  Zn  allow  us  to  make  the  further  computational 
simplification  of  adopting  the  Faber-Ziman  substitutional  model • 5  The 
parameters  used  in  the  computations  and  listed  in  Table  I  are  thus  appropriate 
to  the  a-MgZn  alloys. 

Pseudopotential  matrix  elements  in  the  Born  approximation  and  an  adjusted 
phase  shift  scattering  matrix  element  were  employed.  The  phase  shift  expanded 
form  was  assumed  to  include  only  s  and  p  components  (from  consideration  of  the 
atomic  structure  of  Mg  and  Zn)  and  the  phase  shifts  were  adjusted  to  yield  the 
observed  magnitude  of  p  and  to  satisfy  the  Friedel  sum  rule.  The  Mg  and  Zn 
pseudopotentials  were  taken  from  Harrison. 28  These  two  types  of  matrix 
elements  are  quite  different  in  form.  The  cancellation  near  2kp,  character¬ 
istic  of  pseudopotentials,  does  not  appear  in  the  adjusted  phase  shift 
expanded  matrix  element. 
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The  adjusted  phase  shift  results,  neglecting  saturation  (A  °°)  ,  are  in 
qualitative  agreement  with  the  data,  exhibiting  all  the  observed  features  of 
the  experimental  data  in  a-MgZn  except  for  the  small  minimum  at  about  5K.  The 
phase  shift  results,  including  saturation  with  qpA  =  11.7  (chosen  in 
accordance  with  the  ideas  presented  in  Reference  17),  are  in  remarkable 
agreement  with  the  observed  details  of  the  temperature  dependence  of  the 
electrical  resistivity  in  the  a-MgZn  alloys. 24  This  includes  the  room 
temperature  TCR,  the  shape  and  extent  of  the  (T-T^)3/2  region  to  the  right  of 
the  maximum,  the  magnitude  and  position  of  the  maximum,  the  shape  and  extent 
of  the  quadratic  in  T  region,  and  even  the  position  and  size  of  the  minimum 
near  5K. 

We  also  note  that  negative  TCR’s  occur  at  significantly  higher  2kp/kp 
values  than  predicted  with  the  extreme  backscatter ing  approximation  and  are 
consistent  with  the  data  of  References  24-26.  Furthermore,  the  adjusted  phase 
shift  results,  including  saturation,  and  with  appropriate  values  of2kp/kp,  are 
in  good  qualitative  agreement  with  p(T)  in  the  other  low  resistivity  amorphous 
alloys22 , 23, 25, 26  which  have  been  studied.  On  the  other  hand,  although  the 
Born  approximation  results  are  within  a  factor  of  two  of  the  observed 
magnitude  of  p,  they  fail  to  exhibit  the  temperature  dependence  of  p  observed 
in  the  low  resistivity  amorphous  alloys. 

The  following  conclusions  may  be  drawn  from  the  above  results: 

(a)  The  diffraction  model,  with  appropriate  matrix  elements  and 
incorporated  phonon  ineffectiveness  effects  at  small  qA,  can  explain  the 
observed  temperature  dependence  of  p  in  low  resistivity  (p  <  100  yftcm) 
amorphous  alloys.  Qualitative  agreement  with  experiment  is  obtained  if 
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saturation  effects  are  neglected.  This  latter  is  to  be  contrasted  with  the 
situation  in  high  resistivity  amorphous  metals,  where  not  even  qualitative 
agreement  with  the  observed  temperature  dependence  of  p  can  be  obtained  if 
saturation  effects  are  not  included  in  the  diffraction  model. 

(b)  Saturation  or  Mooij  effects  are  important  even  for  resistivities  as 
low  as  50  yftcra.  Moreover,  the  reduction  in  electron-phonon  interaction 
(phonon  ineffectiveness)  can  be  adequately  represented  by  the  classically 
derived  Pippard  function^  or  even  a  "sharp  cutoff"  form.  Apparently  a 
consistent  procedure  based  upon  a  generalization  of  the  diffraction  model, 
which  assumes  that  saturation  effects  in  elastic  scattering  are  negligible, 
describes  the  temperature  dependence  of  the  resistivity  equally  well  for 
alloys  whose  resistivity  is  50  yftcm  or  150  yftcm.  This  is  particularly 
remarkable  because  it  seems  obvious  that  localization  effects  must  eventually 
lead  to  a  failure  of  this  simple  scattering  picture  and  150  yftcm  corresponds 
to  A  of  the  order  of  interatomic  spacings. 

(c)  There  is  strong  evidence  that  the  Born  approximation  is  not 
appropriate  for  studies  of  the  temperature  dependence  of  p  in  many  disordered 
metals.  Support  fovr  this  is  found  in  the  studies  of  p  in  a-CuSn  by  Frobose 
and  Jackie, 10  the  investigation  of  monovalent  liquid  metals  by  Young,  Meyer, 
and  Kilby, 36  the  study  of  liquid  Ca ,  Sr,  and  Ba  by  Ratti  and  Evans , 38  in 
addition  to  the  present  results. 
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